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Abstract 

We are interested in modeling the Darwinian evolution resulting from the inter- 
play of phenotypic variation and natural selection through ecological interactions, in 
the specific scales of the biological framework of adaptive dynamics. Adaptive dy- 
namics so far has been put on a rigorous footing only for direct competition models 
(Lotka-Volterra models) involving a competition kernel which describes the competi- 
_C tion pressure from one individual to another one. We extend this to a multi-resources 

chemostat model, where the competition between individuals results from the sharing 
of several resources which have their own dynamics. Starting from a stochastic birth 
and death process model, we prove that, when advantageous mutations are rare, the 
population behaves on the mutational time scale as a jump process moving between 
equilibrium states (the polymorphic evolution sequence of the adaptive dynamics lit- 
erature). An essential technical ingredient is the study of the long time behavior of a 
chemostat multi-resources dynamical system. In the small mutational steps limit this 
process in turn gives rise to a differential equation in phenotype space called canonical 
equation of adaptive dynamics. From this canonical equation and still assuming small 
{N) mutation steps, we prove a rigorous characterization of the evolutionary branching 

points. 
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1 Introduction 



Since the first works of J. Monod [22j and Novik and Szilar [23], [23] (see also |25j). 
biologists have developed procedures which allow to maintain a bacterial population at 
a stationary finite size while, at the same time, the bacteria have a positive individual 
growth rate. The procedure is based on a chemostat: bacteria live in a growth container 
of constant volume in which liquid is injected continuously. This liquid contains nutrients 
which are consumed by the bacteria. We assume that the chemostat is well stirred, so 
that the distribution of bacteria and nutrients are spatially uniform. Moreover, since the 
container has a finite volume and fresh liquid enters permanently, an equal amount of 
liquid pours out containing both unconsumed nutrients and bacteria. This pouring out 
helps regulating the bacteria population size. These chemostats are extremely useful to 
make bacteria adapt to the nutrients environment in such a way that they will increase 
their growth rate. This experimental device is for example used (on large scales) in water 
treatment stations. In this paper we study a chemostat model where different nutrients 
arrive continuously and are simultaneously consumed by bacteria which reproduce and 
die in a stochastic time scale. Bacteria are characterized by genetic parameters which are 
inherited during reproduction except when a mutation occurs. These parameters as well as 
the concentrations of each resource, influence the demographics of the bacteria population. 

Usually, chemostat models are essentially deterministic where both nutrients and bac- 
teria population dynamics are described by coupled deterministic continuous process. This 
point of view is based on the fact that reproduction of bacteria happens at the same time 
scale than nutrient dilution. In this work, we consider that the bacteria population is not 
so large that deterministic approximation of its size can be justified. We develop stochastic 
chemostat models based on the previous work of Crump and O' Young |10j (see also [3j), 
where bacteria dynamics is modeled by a birth and death process whose demographic pa- 
rameters depend on the concentration of a unique resource, the later evolving continuously 
in time. This model has been rigorously studied in Collet-Meleard-Martinez [9J 

The literature on multi-resource chemostats is extensive but mostly specialized to very 
precise models. In a more general framework as we consider here, we refer for example to 

[25], m 

Our goal is to study the Darwinian evolution of the genetic parameters of the bacteria 
in the chemostat. We will study how mutations and competition for nutrients lead to their 
progressive adaptation. The study of adaptive dynamics has been developed in the last two 
decades, in particular by Metz, Geritz and coauthors |18[I19]. and Dieckmann and Law [12] . 
This theory emphasizes the connections between ecology and evolution: competition for 
resources strongly influences the selection of advantaged bacteria. The theory was put 
on a rigorous mathematical footing by Champagnat and Meleard and coworkers [5j HI [8] 
with a probabilistic approach (see also [131 [21 [6] for a PDE approach). These results 
are based on the combination of large population and rare mutation scalings, allowing to 
describe evolution as a succession of mutant invasions. In all these works, all the models 
are either Lotka-Volterra models or models of competition for resources assumed to be at a 
quasi-stable equilibrium. This amounts in all cases to assume direct competition between 
individuals. 

Our model is more realistic and takes explicit dynamics for several resources into ac- 
count. This leads to competitive interactions between bacteria, driving selection, with a 
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more complicated nonlinearity than the Lotka-Volterra direct competition (although links 
exist under specific parameter scalings, see [21]). The modeling of adaptation for multi- 
resources chemostat has never been studied from a probabilistic point of view and only 
recently with a PDE's approach (see [20J). Our goal is to study adaptive dynamics in these 
models using the approach of [SJE] combining large populations, rare and small mutations. 
Let us emphasize that adding resources dynamics makes the mathematical analysis more 
complicated and difficult than for the Lotka-Volterra model. 

We first introduce the model in Section [2j We construct a stochastic multi-resource 
chemostat model which couples deterministic and stochastic dynamics. The bacteria dy- 
namics follows a birth and death process with reproduction due to resources consumption 
and mutation. The resource process is deterministic between birth or death events in the 
bacteria population, with coefficients depending on the composition of the population. 
Then we introduce multi-resources chemostat deterministic systems as large size approxi- 
mations of these probabilistic models. 

The core ingredient in our proofs is a stability result on deterministic chemostat sys- 
tems, seeing the stochastic process as a sort of perturbation. In Section [3j we thus study 
the long time behavior of such deterministic nonlinear systems and prove the convergence 
to a unique equilibrium as soon as some non-degeneracy assumptions hold true. The proof 
is based on several Lyapunov functionals which considerably extend the usual Lyapunov 
functional for chemostat systems (see [25J). 

In Section [4j we study the long time stability of the stochastic process, viewed as an 
approximation of the deterministic chemostat system. We prove that all traits with zero 
density at equilibrium actually go extinct after a time of the order the logarithm of the 
population size. We also prove that the time of exit from a neighborhood of the equilibrium 
grows exponentially in the population size. This kind of result follows classically from 
large deviation estimates, which we prove here in the non standard situation of perturbed 
resources dynamics, using the Lyapunov functionals of Section [3] 

We finally study in Section [5] the individual-based process on the evolutionary time 
scale in three steps: First in Subsection 5.1 we consider a large population and rare mu- 
tation scaling but we do not let the size of each mutation go to 0. This ensures that the 
stochastic model on the mutation time scale (evolutionary time scale) converges to a pure 
jump process describing the successive invasions of advantageous mutants. This process, 
called Polymorphic Evolution Sequence (PES), generalizes the TSS introduced in [19J and 
the PES for Lotka-Volterra models whose existence has been proved in [8]. The difficulty 
consists in extending these results to our more complicated model, which is done in Ap- 
pendix [A] Second, we introduce a scaling of small mutations in the PES, letting the size 
of each mutation go to 0. We prove in Section 5.2 that, in this limit, the dynamics of 
co-existing traits is governed by a system of ODEs extending the canonical equation of 
adaptive dynamics (see |12[ [8]). Finally, from this canonical equation and still assuming 
small mutation steps, we prove a rigorous characterization of the evolutionary branching 
points. The only result which needs a different approach from the Lotka-Volterra case [8] 
is the branching criterion, where we use the results of Section [3] to prove that coexistence is 
maintained after evolutionary branching and that the distance between the two branches 
increases. The details of the proof are given in Appendix |Cj after giving useful results on 
the sign of the fitness function (Appendix [Bj , which governs the possibility of invasion of 
a mutant trait in a resident population at equilibrium. 
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2 The Model 



2.1 The stochastic model 

We consider an asexual population and a hereditary phenotypic trait. Each individual i is 
characterized by its phenotypic trait value x % E X, hereafter referred to as its phenotype, or 
simply trait. The trait space X is assumed to be a compact subset of The individual- 
based microscopic model from which we start is a stochastic birth and death process, with 
density-dependence through a reproduction depending on the resources in the chemostat. 
There are r different resources. Resources are injected in the fluid at the constant rate 
1. Their concentrations decrease first because of the linear pouring out of the fluid from 
the chemostat and second by their consumption by the bacteria. We assume that the 
population's size scales with an integer parameter K tending to infinity while the effect 
of the individual resources consumption scales with ^. This allows taking limits in which 

another crucial scale will be the mutation 
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individuals are weighted with In Section 
amplitude a. 

We consider, at any time t > 0, a finite number of individuals, each of them holding 
trait values in X. Let us denote by (x\, . . . ,x n k) the trait values of these individuals. 
The state of the population at time t > 0, rescaled by K, is described by the finite point 
measure on X 

1=1 

where 5 X is the Dirac measure at x. This measure belongs to the set of finite point measures 
on X with mass 1/K 

M K = j^XJ^; n > 0, x u ...,x n e # j . 

Let (u, f) denote the integral of the measurable function / with respect to the measure 
v and Supp(z^) denotes its support. Then (v* , 1) = and for any x € X, the positive 
number (vf, l{x}) i s called the density at time t of trait x. 

The concentrations of the r resources at time t are described by a r-dimensional vector 

n K (t) = «(*),••• ,R?(t))eW + . 

We introduce the following demographic parameters. 

• An individual with trait x reproduces with birth rate given by 



j>k(*)flf, (2.2) 



k=l 

where r?fc(x) > represents the ability of a bacteria with trait x to use resource 
for its reproduction. 
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• Reproduction produces a single offspring. With probability l-fiR-?(4 the newborn 
holds the trait value x of the parent. 

• A mutation occurs with probability fix p(x) with p(x) > and affect the trait of 
the descendant involved in the reproduction event. The trait of the descendant is 
x + h with h chosen according to m a (x, h)dh, where a G (0, 1] scales the mutation 
amplitude. 

For any x the probability measure scales as m u (x J h)dh = ^m(x,^)dh, where 
m(x,h)dh is a probability measure with support {h,x + h G X} and such that 
x + ah G X for all h in the support of m(x, •) and for all a G (0, 1) (this last point 
is for example true when X is a convex subset of Mf). In other words, the support 
of m a has a diameter smaller than crdiam(Af). 

• Small \xk means rare mutations and we assume in the sequel that 

lim K fi K = 0. (2.3) 

K— >OD 

• Each individual with trait x disappears from the chemostat at rate d(x) G (0, +oo) 
either from natural death or pouring out with the liquid. Since the fluid pours out of 
the chemostat at rate 1, one could assume d(x) > 1, although this is not necessary 
for our study. 

• The concentrations of resources in the liquid are solutions of the piecewise determin- 
istic equations: for any k G {1, ■ ■ ■ , r}, 

= 9k -R«-R« (± jf %0*)j = 9k - Rk ~ Rk W?, nth (2-4) 

where > represents the k-th resource injection in the chemostat. The term —Rf£ 
represents the pouring out of the fc-th resource and the term — (u K ,rjk) describes 
the resource consumption by bacteria. Note that any solution to such an equation 
satisfies for any k G {1, ■ ■ ■ , r} the inequality R^(t) < R^(0) A g k - 

The process ((i/f ,R K (t)),t > 0) is a M K x ^-valued Markov process with infinitesimal 
generator defined for any bounded measurable functions <fi from A4 K x (M+) r to R and 
v = if Y!i=i $xi and R = (R u ■ ■ ■ , R r ) by 

L K ^v, R) = | (u + |, i?) - </>(v, R)j (1 - hkp{x)) Vk(x) Rf}j 

J R e i^{ U + ~^K L,R ) ~^^ R ^j VKP(x) (^2vk(x)Rkj m a (x,h)dh 
v-jt,Rj -<f>(v,R)^j d(x)^Ku(dx) 

+ j2ljf^ R ^9k-Rk-RkW,Vk))- (2.5) 
k=i k 



+ 



+ 
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The first term describes the births without mutation, the second term the births with 
mutation, the third term the deaths and the fourth term the dynamics of resources. 



2.2 An example 

We consider a case with two resources and one-dimensional traits having opposite effects 
on the two resources consumption. More precisely let us define the following parameters: 

X = [—1,1], fj,jcp(x) = p, m(x,h)dh = M(0,a 2 ) (conditioned on x + h G X), r = 2 (2 
resources), g\ = gi = 1, d(x) = | + ax 2 with a > (minimum at 0), r]i(x) = (x — l) 2 , 
Vi{ x ) = (x + l) 2 . The parameter a measures the impact of trait on mortality. High a 
means strong effect of traits away from on mortality. 



Simulations are given in Figure 2.2 where initially K individuals have the same trait —0.5, 
for three different mortality parameters, a = 1/4, a = 1/2 and a = 1. The upper panels 
represent the time evolution of the density of traits in the population, and the lower panels 
the resources concentrations. We observe two different behaviors. In the four pictures, the 
support of the population process first approaches 0, where mortality is the smallest, but 
next, for a = 1/2 and a = 1, the support of the population process stays close to for 
a long time, while for the two simulations with a = 1/4, the population divides into two 
subpopulations with distinct trait values, but still interacting for the same resources. This 
phenomenon is known as evolutionary branching. 



The two different simulations of evolutionary branching for a = 1/4 (Fig. 2.2 (a) and (b) ) 
illustrate two slightly different ways of branching. In the first one, evolutionary branching 
occurs when the traits in the population are distributed around a value not exactly equal 
to 0, contrary to the second simulation. As will appear in this paper, this explains why, in 
the first picture, the two branches are of very different size at the beginning of evolutionary 
branching, while they are of similar size in the second picture. 



Fig. 2.2 (c) does not show evolutionary branching, but the width of the trait distribution 
around is wider than in the fourth picture, where a is bigger. As will appear below, the 
third picture actually corresponds to a critical case with respect to evolutionary branching. 



Finally, in Fig. 2.2 (d), we also observe an interesting situation in the first phase of evolu- 
tion, where the trait distribution in the populations is approaching to 0. There is a jump 
in the support of the trait distribution around time t = 300, which can be understood as 
follows: at this time, a mutant^] just appeared by chance, with a much bigger trait than 
the largest trait in the population just before. This mutant survived and produced a large 
descendence with a significantly smaller death rate than the rest of the population, which 
was therefore competitively disadvantaged, and went extinct fast due to this competitive 
pressure. 

Our goal in this paper is to give a description of these pictures under a specific scaling of 
the parameters of the individual-based model, and to give some conditions under which 
evolutionary branching appears. The evolution of the population results from an instanta- 
neous trade-off between death rate minimization and the birth rate maximization (through 
a better consumption of resources) of the individuals. 



1 Actually two different mutants traits, if one looks closely at the simulation. 
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1000 2000 3000 4000 5000 




1000 2000 3000 4000 



1000 2000 3000 4000 5000 



1000 2000 3000 4000 



(a) K = 300, p = 0.1, a = 0.01, a = 1/4 (b) K = 300, p = 0.1, a = 0.01, a = 1/4 





1000 2000 3000 4000 



(c) K = 300, p = 0.1, a = 0.01, a = 1/2 (d) K = 300, p = 0.1, a = 0.01, a = 1 



Figure 2.1: Simulations of the individual-based model for three different values of the 
parameter a. Upper panels: time evolution of the trait density in the population. Lower 
panels: time evolution of the resources concentrations. 
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2.3 The process at the ecological time scale 

The ecological time scale is the birth and death time scale in which the process is defined. 
The evolution takes place on the longer time scale of mutations. 



The following properties will be assumed in the sequel. 

Lim.tr.hitz continuous on X. 

(2.6) 



The functions rji, . . . ,rj r and d are C 2 on X and p is Lipschitz continuous on X. 



Since X is compact, these functions are lower bounded by positive constants. 

The function m(x, h)(and thus m cr (x, h)) is Lipschitz continuous on X x Mr. 

In addition, for any a G (0, 1], there exists a function m a : M — > M + such that ^ 7) 

m a (x, h) < rh a {h) for any x € X and liGK and I m a (h)dh < +00. 



The initial conditions satisfy supE((z/^, 1)) < 00 ; sup i?f"(0) < 00. (2. 

K K,k£{l,...,r} 



For fixed K, under (2.6)-(2.7)-(2.8), the existence and uniqueness in law of a process 
on D(R_)_,.M x with infinitesimal generator L K can be adapted from the one 

in Fournier-Meleard [15] or |4]. The process is constructed as solution of a stochastic 
differential equation driven by point Poisson measures describing each jump event plus a 



drift term describing the resource dynamics and Assumptions (2.6) and (2.7) prevent the 



population from exploding since resources concentrations are bounded. 

Let us recall the construction. Let Ni(ds,di,d0), N2(ds,di,d0,dh) and N^(ds,di,d9) be 
independent Poisson point measures on M + x N x IR + , IR + x N x R + x R e and M+ x N x R+ 
respectively, and with intensity measures qi(ds,di,d9), q2(ds,di,d0,dh) and q^(ds,di,d8) 
respectively, where 

qi(ds,di,d9) = q 3 (ds,di,9) = l{ s > ,e>o}ds ^2^ k (di)j d8 

and 

/ 00 \ 

q2(ds,di,d6,dh) = l^ s>0 ^ >0 jds ' s ^8k{di) d9fh a {h)dh. 

\k=i ) 

For all v £ M K , we define Xi{v) £ X for all i > 1 such that v = (^«>i ^0)) an ^ 

xi(v) z< x 2 {v) z< ■■■ r< x n{v) (v), 

where n(v) := K(i>, 1) and ^ is an arbitrary total order on X C (e.g. the lexicographic 
order. 
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Now, let us consider the equation 



.K , K 



v n + 



JN JO 



+ 



JN JO 

t 



JI Jo 



^i(^) ]1 -{i<n(^), 8<d(x t (vK_))} N 3(ds,di,d8) 



(2.9) 



coupled with the equations (2.4) for resources with Rk(0) > 0. 



This system of equations admits a unique solution, up to a time of accumulation of jumps. 
The next proposition shows that this time is infinite and gives uniform moment estimates. 



Proposition 2.1 Assume (2.6) -(2. 7) -(2.8), then the process does not explode in finite 
time almost surely and for all t > 0, 



sup sup E 

K>1 t>0 



fc=l 



< oo. 



(2.10) 



Proof Since < Rf(t) < g k V ||-Rf (0)||oo =: R* k for all t > 0, the individual birth 
rate in the population is always smaller than Ylk=i ^fc^fc' wnere Vk := ll^fclloo- Therefore, 
for all t < 0, K(u^ , 1) is stochastically dominated by a pure birth (Yule) process in Z + 
with transition rate k Ylk=i Vk^-k f rom k to k + 1 (this process can be easily explicitly 
constructed from the point processes N±, N2 and N3). The Yule process is a.s. finite for 
any t, so that the process is well defined for all times. In addition, we get that 



k=\ k=\ Jo Jx \k=i 

+ E f (9k - Rkis) - Rg{8){v?, Vk }) ds + M t , 
,._ n Jo 



d(x) 1 v s (dx)ds 



with 



Mt 



1 

K 
1 

K 
1 

K 



Jn Jo 

t 



1 



_))R«(s)}Ni(ds,di,d6) 



-{i<n{vf_), 0<m( a: 4 I /JL),ft)^p(a: i (^))EI = i%(^(^il))KfW}^ 2 ^ S '^' d6 '' d/t ) 



1 {i<n(^i), 0<d(xi(uf_))} N s(ds,di,d9), 



where TV, = iVj — i = 1,2,3, are the compensated Poisson processes. The process 
(Mt, t > 0) is a martingale and thus 



E 



fc=l 



E 



<^,l)+5>f(0) 



k=l 
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(2.10) then follows from the fact that d(-) > d > by Assumption (2.6) and Gronwall's 
lemma. In addition, because of the shift invariance of the Poisson point measures N\, N2 
and N3, the process (u^ ,H K (t))t>o is strong Markov. □ 



The next result shows that mutations cannot occur on bounded time intervals, since the 



mutation time scale t/K/j.^ tends to infinity by (2.3). 

We define T^ nt as the first mutation time of the population process (f/^, t > 0). 



Corollary 2.2 Assume (2.6)-(2.7)-(2.8). Then for allrj > 0, there exists e > such that 
for allt>0 

t t + e 



limsupP ( a mutation occurs on 

K— >oo 



K^ K ' KjJL K 



< 77. 



Proof The proof is based on a coupling inspired by the previous proposition. By the 
Markov property, it is sufficient to prove the result for t = 0. Let vt be the solution of 



V t = VQ + 



1 

K 



t p POO 



10 jnjo 
coupled with the equation 



() Jr /n 5 ^s-) 1 {i<n(u s _), 0<(l-^p(a :j (P s -)))Efc=l%(^(^-))«fc(s)} iVl ^ S, ^' d0 ) 

^ 4 (i/ s _)l{i<n(P s _), e<d(x x (i> s ^))} N 3(ds, di, dO) 



dR k (t) 
dt 



g k - R k - R k (i>s,ilk) 



Comparing with (2.9) we see that vt = vt an d R k (t) = Rk{t) for all t < T^ ut . 
Now, let us define the Z + -valued process 
i*t r poo p 

II I !{i<n(i5 s _), 8<mv(h)ii K rrig} N 2(ds,di,d9,dh), 

Jn Jo Jr 1 

where fj = maxjsup^g^ rjk(x); 1 < k < r} and g = max{||i?fc(0)||Loo V gk', 1 < k <r}. Then 
^rnut — Ta, where Ta is the first jump time of (At,t > 0). Since iV~2 is independent of 
(z/, R), we have 



T K 

mut 



Kuk 



> P [T A > 



K 



E 



e/Kfif 



exp 



> 1 - K/j, K rfjg I E{uf, 1) ds 



n(uf)ii K rr]gds 



and Lemma 2.2 follows from the boundedness of the moments (cf. Proposition 2.1). □ 



10 



2.4 Convergence to deterministic chemostat systems when K — > +00 

In this section we study the large population and rare mutation approximation of the 
process described above when the initial measure has the finite support {x±, • • • , x n }. The 
limit is deterministic and continuous and the mutation events disappear. The next result 
is a simple but useful first step to characterize the dynamics between mutation events. 

We introduce the following chemostat (coupled) system, denoted by CH(ra, zi, ■ ■ • , x n ), 
solved by (u(t),R(t)) € R^ +r : For % = 1 • • • n, for k = 1 • • • r, 



d ( x i) + ^2m(xi) R k ); 



m = it. 

fc=i 

Rk = 9k - Rk 1 + ^ Vk(xj) Uj 



(2.11) 



Theorem 2.3 Let x\, . . . ,x n be distinct points in X . Assume that Uq — Y2i=i u o ' ^hh 
such that Uq' 1 — > u l and Rt?(0) — > Rkfi in probability, where Uq and Rkfl are deterministic, 
nonnegative numbers. Then, for all T > 0, 



lim P(T* ttt <r)=0, 



and 



sup [^\(^ K a Xi }-Ui(t)\ + J2\ R k^- R ^\) ^° 



0<t<T 



k=l 



in probability as K — > +00, where (ui(i), . . . , u n (t), Ri(t), . . . ,R r (t)) is the solution of the 
chemostat system CH(n, x%, . . . , x n ) with initial condition itj(0) = u l and i?fe(0) = Rk,o- 

Proof Since Kfix 0, the fact that linift'_>. +00 V(T^ ut < T) = follows trivially 
from Corollary |2.2| Therefore, in order to prove the second part of the result, it suffices 
to prove it for the population dynamics obtained by setting to zero the birth rates with 
mutation. In this case, the model reduces to a birth and death Markov chain in for 
{{v^ , la*) j 1 ^ i < n )t>o, coupled with piecewise deterministic dynamics for the resources. 
Then the second part of the result can be proved using standard techniques from [14, 
Ch. 11]. The only difficulty comes from the fact that the birth and death rates and the 
vector fields of the resources dynamics are only locally Lipschitz functions of the state of 
the process. Since the limit function (ui(t), . . . , u n (t), R±(t), . . . , R r (t)) takes values in a 
compact set, the difficulty can be easily solved by regularizing the transition rates out of 
a sufficiently large compact set. □ 
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3 Asymptotic behavior of the deterministic chemostat sys- 
tem 



3.1 Assumptions and statement of the results 



In order to study the long time behavior of the system (2.11), we need some additional 
assumptions. 



For any x G X, Vk( x )9k > d(z 



(3.1) 



k=l 



For any n > 1 and any distinct x\, ■ ■ ■ , x n € X, Equation (3.3) has at most one 
solution (ui, • • • , u n ) £ W±., where 



o 



1 < i < n. 



(3.2) 
(3.3) 



Assumption (3.1 ) means that when resources are maximal the population process is super- 



critical. Thus, (3.1 ) prevents the individual-based model to become extinct too fast. It also 



ensures that the trivial equilibrium (0, • • • , 0, gi, ■ ■ ■ , g n ) of the deterministic system (2.11 ) 
is unstable. 

Since the equilibria (u, R) of the chemostat system CH(n, x±, . . . , x n ) are given by canceling 



the right hand side of (2.11 ), they must satisfy 

Rk 



9k 



1 + Y,]=iVk(x j )u j 
and for all i either ui = or d{xi) = Ylk=i Vk( x i)Rk- 



Therefore Assumption (3.2) implies that, for all / C {l,2,...,n} there is at most one 
equilibrium (u,R) of (2.11) such that Ui = for all % £ I and iZ, > for all i E I. In 



particular, we will make use of the following consequence of (3.2 ): if (u, R) is an equilibrium 



of (2.11 ) and v is a vector of such that V{ = implies m = for all 1 < i < n and for 
any k 

n n 

Y^m&jhj = ^2vk{x j )v j , 

then v = u. 



Proposition 3.1 (i) Assumption (3.2 ) is satisfied as soon as: For all distinct x%, ■ ■ ■ , x r +i £ 
X , the vectors 

( m( x i) \ ( Vr{xi) \ ( d(xi) \ 

: : ,| : (3.4) 

\qi(x r+1 )J \r] r (x r+1 )J \d(xr+i)/ 
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are linearly independent, and for all distinct xi, . . . , x r , the vectors 



(3.5) 



\r]r(Xr) / 



are also linearly independent. 



(ii) Conversely, Assumption (3.2) implies that, for all distinct x\,...,x n such that (3.3) 
admits a solution with all positive coordinates, the vectors 



Vect 



'rji(xiy 



( T]r(xi)\ 



\r]r(x n )J 



(3.6) 



are linearly independent. In particular, (3.3) admits no solution in (0, +oo) n ifn>r. 



Proof Let us first assume that n > r + 1 and fix x\, . . . ,x n distinct. In view of (3.4), 
the system of equations 



d{xi) - ^2 Vk(xi)R k = ; 1 < i < 



n 



k=i 



has no solution (R\, . . . , R r ). Hence the system (3.3) has no solution. 

If n < r, consider two solutions u and v! of (3.3) and define the vector R = [R\ 

9k 



by 



Rk 



1 + Si Vk{Xi)Ui 
and the vector R' similarly in function of v! . Then 

\ 2 



VI < k < r, 



R, 
(3.7) 



RkR'k 

9k 



Yl J2( u i ~ ^-^-^Vkix^^Vkixi)^ - u'i) 

k j 9k i 

I ~ u i Vk(xj)Rk - — YV- Vk(xj)R' k I V r) k {xi)(ui - u-) 
k \ 9k j 9k j 



Y( u i - u'i)J2m(xi)(Rk - R' k ), 



where we have used (3.7) for (R, u) and (R' , u'). Since rj k (xi)R k = ^2 k Vk(xi)R' k = d(xi 
for all i, the previous quantity is 0. Therefore, 

y^Vk(.Xi)(uj - u'j) = 0, Vfc G {!,..., r}. 
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Since n < r, the vectors 



fm(xi)\ A/i On) 



W^l), 



are linearly independent under condition (3.5), which ends the proof of (3.2) by implying 
that Ui = u[, Vi. 

Assuming that Point (ii) does not hold, and letting x±,...,x n be distinct traits and 
U\, . . . , u n be a solution to (3.3) such that 



atirjk(xi) + . . . + a n r) k (x n ) = 
for all 1 < k < r and for some (otx, • • • , Q n ) 7^ 0. Then, for e close enough to 0, the vector 



(v\, . . . , v n ) belongs to (0, 00 ) n and is another solution of (3.3 ), where V{ = Ui + □ 



The two assumptions (3.1)-(3.2) imply our main result on the large time behavior of the 
chemostat systems of ODEs. 

Theorem 3.2 Assume (3.1) and (3.2). For all n > 1 and all distinct x±,--- ,x n E X, 



there exists (u,R) in (M+) n+r such that any solution (u(t),R(t)) of the system (2.11) with 
U{(0) > for any 1 < % < n, converges to (u,R). In addition, (u,R) is the unique 



equilibrium of the system (2.11) satisfying, for all i such that Ui = 0, the (weak) stability 
condition 



d{xi) + Vk(xi) R k < 0. 



(3.8) 



k=l 



Before giving the proof, we define two notions of great importance in the sequel. Writing 
x = (xi, . . . , x n ), we denote by 

u(x) = (ui(xi, . . .,x n ), . . .,u n (xi, . . .,x n )) 

the vector of equilibrium densities of the previous theorem and by 

R(x) = (Ri(xu . . .,x n ), . . .,Rr(xu . . .,x n )) 

the corresponding equilibrium resources concentrations. 

Definition 3.3 We say that the traits x±, ■ ■ ■ ,x n coexist if the quantities ?2i(x), . . . , n n (x) 
are all positive, where x = (x\, . . . ,x n ). For all n > 1, we denote by D n the domain of 
coexistence of n traits: 

T^n = {(xi, ■ ■ ■ ,x n ) G X n : x%, . .. ,x n coexist}. 

Note that T>i = X. 

For distinct traits xi, . . . ,x n , we also define the invasion fitness of a new trait y as the 
function 



f(y,xi, 



in) = -d{y) + ^2 Vkiv) R k 



(3.9) 



k=l 
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Remark 3.4 We can interpret Assumption (3.2) and Proposition 3.1 with this vocabu 



lary: (3.2) implies that there exists at most a single vector of population densities where 



xi, . . . , x n coexist for all distinct x\, . . . , x n 6 X, and Proposition 3.1 (ii) means that when 
xi, . . . ,x n coexist, we must have n < r. In particular, T> n = if n > r. 

If the trait y is interpreted as a mutant trait trying to invade the resident populations 
of coexisting traits x\, . . . ,x n , the terminology of fitness refers to the possibility of in- 
vasion of the mutant trait. Indeed, it follows from Theorem 3.2 that u(x\, . . . ,x n ,y) = 



(«i(x), . . . , n n (x), 0) iff f(y; x\, . . . , x n ) < 0, and that if f(y; x\, . . . , x n ) > 0, the last equi- 
librium is locally unstable, since the Jacobian matrix of the system at this point obviously 
admits f(y; genvalue. 



3.2 Proof of Theorem [3721 

The proof of this results was sketched in [7]. We shall give here a complete and detailed 



proof. The idea is to construct a Lyapunov functional for the system (2.11 ) in three steps. 



Step 1. Lyapunov functional and stable equilibrium for a reduced system: We 



consider the quasi-stable approximation of the system (2.11) obtained by putting at each 



time t the resources at the equilibrium associated with the population densities at time t. 



Ui 



v rik(xi)g k 



Vi G {1, . . . , n}. 



This system admits the Lyapunov functional 



(3.10) 



F(ui, ...,u n ) = s ^d{xi)ui - ^ g k log 1 + ^^rj k (x j )u j . 



(3.11) 



One easily checks that 



d 
dt 



This is a strict Lyapunov functional since dF(u(t))/dt < 0, and dF{u{t))/dt = if and 



only if u{t) is an equilibrium of the system (3.10). Clearly, 



F(u) — > +oo when |u| —> +oo with u G , 
since d(x) > for all x € X. Therefore, for any initial condition u(0), the solution u(t) 



of (3.10) converges to an equilibrium when t — > +oo. In addition, F is a convex function 
as 

d 2 F v-^ 9k t \ t \ 



duiduj 



+ T,iVk(xi)uiY 



Hence 



dUiduj ^(i + E 



9k 



J^Vkixi)^) >0, (3.12) 
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thus the convexity. In general F is not strictly convex as (3.12 ) can vanish for nonzero £ if 
n> r. However we can still prove that there is a unique local minimum. Indeed, consider 
any critical point u of F. Define I the set of indices i s.t. m > 0. Then for any i E / one 
has that 



d(xi) 







We denote R 1 the subspace of W 1 defined by £ E M / iff & = for any i £ I. At the point 
u, F is strictly convex in M 7 . Indeed if not, one could find £ E R s.t. ??fc(xj)£* = for 
any fc. 

On the other hand for e small enough u + e£ also belongs to Wt and one also has that 



d{xi) - X) i 



9k f]k{Xi) 



0. 



which would violate the Assumptions (2.6) and(3.2). This shows that if u is a local 
minimum of F then no other local minima may exist on M 1 , as F is convex over R and 
strictly convex at u. 

Since F(u) —> +oo when \u\ —> +oo with u E M™ , F has at least one global minimum. 
Choose a global minimum u such that |/| is the largest and assume that there exists 
another one u', defining another set I'. Since F is strictly convex in Mr, one necessarily 
has that I' <f_ I. In addition by the convexity of F any u 9 = 9u + (1 — 9)u' is also a global 
minimum. However for 6 E (0, 1), one has that uf > for any i E / U I' which is strictly 
larger than I. Hence we obtain a contradiction and F has a unique global minimum. 



Therefore, F admits a unique global minimizer in the closed, convex set 
Let us denote by R the vector with coordinates 



1, denoted u. 



Rk 



9k 



1 + J2i=lVk(Xi)Ui 



Kk<r. 



Then, the vector (u,R) is an equilibrium of (2.11). Since u is a glo bal minimum of F on 
M™, for any i such that = 0, one must have JJjj- > 0. This yields (3.8). 



Let us check that (u,R) is the only equilibrium of (2.11) satisfying this property. This is 
equivalent to checking that u is the only equilibrium of (3.10) such that 



k=i ^ 



9k rjk{Xi) 



jVk(xj)uj 



< 



for all i such that U{ = 0. Since u is an equilibrium of (3.10), note also that 

9kr]k(xi) 



-d(xi) + ^2— ^ 
k=i z -'- 







for all i such that m > 0. 
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Let us consider two such equilibria, u 1 ^ u 2 . Then, adapting the computation for the 
convexity of F 

o > w [-dM + v - 9 " m{x ;\ 2 ) 

T V * i + Ej^^iKv 

gkrik(xi) \ 



E/-1 -2\ | 7/ \ i 9kVk( x i) \ 

( ( " , -"' ) r (i,)+ ?i + E J -^K 2 J 

i V fc A 



3fc ??fe^ij 



Eft 



.9/,: 



ul-u}) 



This cannot hold unless 



X)ffc(xi)(u?-«J) = 0, Vl<A:<r. 



But this would imply that i? 1 = i? 2 (with obvious notations). Defining J := {i : 
^2kVk(xi)Rl = d(xi)}, since u 1 and u 2 are equilibria of (3.10), we would then have 



for all i £ J, and we would obtain a contradiction with Assumption (3.2) 



applied to (iCi)jgj. Hence (u,R) is the only equilibrium of (2.11) satisfying (3.8) 



Step 2. A degenerate Lyapunov functional for the system (2.11): Let us fix 



(u(0),i2(0)) G (0,+oo) n+r and consider the solution (u(t),R(t)) of ( |2.11) . 
We define 

n r 

. . . ,u n ,R!, . . .,R r ) = }Xuj - UilogUi) + ^(Rk - Rk^ogR k ). 



i=l 



k=l 



Then, for any solution of ( |2.11| ) with (n(0), i?(0)) G (0, +oo) n+r , one has 



£(«(*),#(*)) = £( 



- Ui 



-d(xi) + ^2vk(xi)Rk 



(3.13) 
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or 



dt 



G(u(t),R(t)) = - Ui)^2rik(xi)(R k - R k ) + E( 



[Ui — Ui, 
% k 



Ui - Ui 



-d(xi) + y^ j rj k (x i )Rk 



rik(xi)ui) 

k k \ i 

- ^(Rk - Rk)y^r)k( x i)( u i ~ ^i). 



(3.14) 



Equation (3.8) implies that the second term in the r.h.s. is non-positive. Therefore, 



d G(u(t),R(t)) < ^2 Rk Rk Rk (Rk ~ Rk) ( 1 + E ^)u t 



dt 



+ E ^ R^ ( 9k ~ Rk ^ + E Vk(xi)Ui) 



2 7 ^2 I 5 



(3.15) 



since the second term of the r.h.s. is zero. Therefore, G is a Lyapunov functional for the 



system (2.11), degenerate in the sense that, in view of (3.14), its derivative could vanish 



when R k = Rk but U{ ^ Ui for some i. 

Note that G is convex (strictly in the Rk and in the Ui for which U{ 7^ 0) and G(u, R) — > +oo 
when \u\ + \R\ — > +oo in M" +r . As a consequence, the function u(t) is bounded, say by U. 

Since Rk < 9k — Rk, Rk(t) is uniformly bounded in time, and by Lyapunov's Theorem, 
R(t) — > R when t — > +oo and thus 



Now, let us define 



a := inf inf Ru(t) > 0. 
*>o k y ' 



!:={%■ ^2vk(xi)Rk = d(xi)}. 



Note that if Ui > then i £ I but that / could be larger. Since (u,R) satisfies (3.8), 
-d(xi) + Yjk^kixijRk < for all i /, i.e. 



P := inf d(xi) - Y] r)k(xi) Rk 
k 



> 0. 



Then, using (3.14), one can find a more precise expression than (3.15) 
d -G{u{t),R{t)) < -p 



dt 



G(u(t),R(t)) <-/?I>-E R * k? ( 1 + E W 



(3.16) 
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Step 3. Second Lyapunov functional for the system ( |2.11[ ): We can partly correct 
the problem that G is a "degenerate" Lyapunov functional by slightly modifying it: let 7 
be a small positive number and define 



where 



Then 



H(u, R) = G(u, R) + jG(u, R), 
G(u,R) = ^2(Rk ~ Rk)^2m(xi)(ui - Hi). 

k i 

^G(u(t),R(t)) = J2 \9k-Rk(l + ^2^(x J )u J )) Y,Vk(xi){ 

k \ j 
- ^Rki ^2vk(Xi){Ui - Ui 



(3.17) 



Ui - Ui 



+ ^2( R k - Rk)^2vk(xi)ui I -d(xi) + ^2r]i(xi)Ri 

k i V I 

+ Ui r) k (xj)(R k - R k )j . 



(3.18) 



Since (u,R) is an equilibrium of (2.11), the first term in the r.h.s. is equal to 



^ 'R k -Rk) 2 
2R?R k 2 



It follows from the definition of the set / that the third term of the r.h.s. of (3.18) is equal 
to 

^(Rk ~ Rk)^2m(xi)ui -d(xi) +y^ j rji(x i )Ri <C^Ui, 
k i0 V I J igl 

where the constant C depends only on uniform upper bound for the functions Rk(t), f]k{x) 
and d{x). Finally, using Cauchy-Schwartz inequality and the boundedness of u, the last 



term of the r.h.s. of (3.18) is bounded by 



(R k - R k f 
R k 



where the constant C depends only on U and uniform upper bounds for Rk(t) and rj k (x). 
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Combining all these inequalities with (3.16), we obtain 
-H(u(t).R(t)) < - | l-T^S 4 -7C) ^ 



|ff(»(<), < - (l - " 7C) £ 



2 

7 
2 



A: \ i / m 



Ui 





(3.19) 

if 7 is small enough. Lyapunov's Theorem then entails that the set of accumulation points 
of (u(t),R(t)) when t — > +oo is contained in the sub-manifold M of M™ +r defined as the 
set of (vi, . . . ,v n , Si, . . . , S r ) € M™ +r satisfying 

S k = R k , yi<k<r 

n 

y^Vk(?i)(vi-ik) = 0; VI < fc < r, 
i=l 

Wi = 0, Vi J. 

Since iZj = for all i I, this sub-manifold M contains the point (u,R). Since (u,R) is 
an equilibrium, the second system of equations above is equivalent to 

9k - 9k -R k , Vl<k< r. 



l + 12iVk(xi)vi l + ^2iVk(xi)ui 



Therefore Assumption (3.2) applied to the vector of traits (xj)j g / implies that V{ = Ui 
for all i £ I. Hence M is reduced to the point (u, R) and the proof of Theorem 3.2 is 
completed. 

3.3 Some examples 

3.3.1 The monomorphic case 



In case n = 1, we consider a unique trait x. The equilibrium u(x) defined in Theorem 3.2 
satisfies 



1 + r] k (x)u(x) 



Remark that the left hand side is a decreasing function of u(x), and so there is a unique 
solution to this equation. Therefore by the Implicit Function Theorem and Assumption 



(2.6), the function u(x) is a C 2 -function. 

thus g 

Rk(x) 



The resources at equilibrium are thus given by 

9k 



1 +Vk(x)u(x)' 
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3.3.2 The dimorphic case 

In case where n = 2, we consider two distinct traits x\ and x 2 . Then the fitness function 
defined in (3.9) is given by 

Vk{.x 2 )gk 



f(x 2 ; Xl ) = -d(x 2 ) + £ ' ,kK : 2 \ y * ( y (3.21) 
f— ' 1 + mixi u[x-[ ) 



- +Vk(Xl)u(x 1 ) 

It is immediate to observe that / is C 2 with respect to {x\, x 2 ). 

The system CH(2,xi,X2) has two obvious equilibria: E\ = (u{xi), 0, R(x{)) and E 2 = 
(0, u(x 2 ), R(x 2 )). In view of (3.8), E\ is the equilibrium given by Theorem 3.2 if and only 
if f(x 2 ;xi) < and a similar condition for E 2 (and we cannot have both f(x 2 ; X\) < and 
f(xi;x 2 ) < 0). If we have f(x 2 ;x\) > and f(xi;x 2 ) > 0, the equilibrium of Theorem 



3.2 must have positive density coordinates (ui(x\, x 2 ), u 2 (x\, x 2 )). Therefore, we have the 



following result. 

Proposition 3.5 For any x\ ^ x 2 in X , the traits x\ and x 2 coexist if and only if 
f{x\; x 2 ) > and f(x 2 ; x\) > 0. 

Moreover the vector (ui(x\, x 2 ) , u 2 (x\, x 2 )) is solution of the system 

^ 1 + T] k [x\)u\ + rj k (x 2 )u 2 

d(x 2 ) = ±- /frfo)** . (3.23) 

^ 1 + 7? fc (xi)ui + rj k (x 2 )u 2 

3.3.3 The n-morphic case 

When the population is composed of n distinct traits xi, . . . ,x n , the case of coexistence 
can be characterized by induction over n by negating the stability condition (3.8) of Theo- 
rem 3.2 for each trivial equilibrium of the system CH(n, x\, . . . , x n ). By trivial equilibrium, 



we mean any equilibrium with at least one null coordinate for population densities. 
In the case where x = (x\, . . . , x n ) coexist, the equilibrium densities satisfy 



d^) = E i+^ kVk{ 7\- ( V V ? e {!,..., n}, (3.24) 



and the resources concentrations at equilibrium are given by 

* w= i + gJU)M«) ' v * e{1 " -'■>' (3 - 25) 

We have the following regularity result for the function u on T> n . 

Lemma 3.6 For all n > 1, the function u is C 2 and bounded on T> n . In addition, T> n is 
an open subset of X n . 
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Proof Let us first prove that u is C 2 on T> n . For x = (x\, . . . ,x n ) G D n , the equilib- 



rium population densities Uj(x) are the (unique, by Assumption (3.2)) solutions of (3.24). 
Defining the function F = (Fx, . . . ,F n ) from X n x W 1 to R n by 



" i\XXi ■ ■ • > X n , U\, . . . , U n ) — U[Xi) 2_ j 

k=i 1 + ^i=i 



9kVk (xi) 

rjk(xj)ui 



Vi G {1, . . . ,n}, 



we see that u(x) is characterized by the equation F(x, u(x)) = 0. The Jacobian matrix of 
F with respect to (ux, • • • , u n ) at (x, u) = (xx, ■ ■ ■ , x n , ux, ■ ■ ■ , u n ) is given by 



gkVk{xi)f] k {xj) 



J u F(x,u)= Y, 

\k=i (l + Y. j Vk(xj)u j 



l<i,j<n 



For all vector v G 



v*J u F(x, u)v = Y — x2 

k=1 i 1 + J2jVk(xj)uj 



Y Vir lk( x i) 

.1=1 



Hence the matrix J u F(x, u) is invertible if and only if the orthogonal vector space of the 
family of vectors {rjx (x), . . . ,ry r (x)} in W l is {0}, where %(x) = (rjk(xx), ■ ■ ■ ,rjk(x n )). This 
is equivalent to the fact that the vector space spanned by f?i(x), . . . ,rj r (x) is M. n . This 



is implied by (3.6), which is itself a consequence of Assumption (3.2) by Proposition 3.1 
Hence the fact that u is C 2 on T> n follows from the Implicit Functions Theorem, as well as 
the fact that T> n is open. 

Let us now prove that u is bounded on T> n . This property for all n > 1 is equivalent to the 
fact that the function u is bounded on X n for all n > 1. We shall prove this last point. 



In the proof of Theorem 3.2 the equilibrium u(x) is characterized as the unique global 
minimizer on W}_ of the function 



(ux, ...,u n )^ F(xx, . . .,x n ,ux, 



n 

E 

i=i 



d{xi)Ui - Y 9k log { 1 + Y Vk(xj)uj 

3=1 



k=l 



Now F(xx, . . . , x n , 0, . . . , 0) = for all x G X n and 

F(xx, ...,x n ,ux,...,u n )> d{ui + . . . + u n ) - rglog(l + fj(ux + ... + u n )), 
where d := min x( z X d(x) > 0, g := sup 1<fc<r 5' fc and fj := sup k \\r] k \ 



loo. Therefore, the 

function u(x) is uniformly bounded on X n , and Lemma 3.6 is proved. □ 



The coexistence criterion of Proposition 3.5 extends as follows to the n-morphic case. 



Proposition 3.7 For any x G T> n and y G X such that for some i G {1, . . . , n}, 

(pcx, • • • , Xi—XiV-, Xi+l, ■ ■ ■ , x n ) G T> n (this holds for example if y is close enough to Xi). Then 
(x,y) = (xx, . . ■ ,x n ,y) G V n+1 if and only if /(y; x) > and /(a;,; xi, Xi-x, V, ■■■ ,x n 
0. 
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Proof In the first step of the proof of Theorem 3.2 u(x, y) was characterized as the 
unique global minimizer of the functional F defined in (3.11). By Theorem 3.2, we also 
know that the unique global minimizer of F on X n x {0} is u* := (u(x), 0) and the unique 
global minimizer of F on X l ~ l x {0} x X n ~ i+1 is 

u 2 : = (^l( x *)> • • • ,Ui-l(x*),0,U i+1 (x*), . . . ,u n (x*),Ui(x*)\ 

where x* = (x\, . . . , Sj_ i, y, • ■ ■ , x n ). Now, since the derivative of F with respect to 
the n + 1-th coordinate at u* is 

-/(y;x) <0 

and the derivative of F with respect to the i-th coordinate at u?j is 

the unique global minimizer of F must have positive coordinates, which is the definition 
of (x,y) 6 V n+1 . □ 



4 Long time stability of the stochastic process 



Theorem 3.2 gives the asymptotic behavior of the deterministic chemostat systems. When 
K is large, one can see the individual-based model as a stochastic perturbation of the 
deterministic chemostat system. One may then wonder to what extent this perturbation 
modifies the long time stability of the system. The next result answers this question and 
is going to be useful in Section [5] 

Let x\, . . . , x n be distinct points in X and e > and denote by M. the set of finite measures 
on X. We define 



B e (x) = |(i/,R) eMxW : Supply) C {xi,...,x n }, V* \{v, l {xi} )-Hi(x)\ < e, ||R-R(x)|| < 

Theorem 4.1 Let x\, . . . , x n be distinct points in X. Assume that Vq = X^ILi u o' l $xi such 
that (uq , H K (0)) a.s. belongs for K large enough to a compact subset S of (0, +oo) n x W . 

(i) For all e > 0, there exist T £ , V e > such that 

lim IPfvt G [T e , e VeK A T^ ut ], (yf , R K (t)) € B e (x)) = 1. 

A— >oo \ J 

(ii) Define I = {i : Uj(x) = 0} and T^, t = inf{i > : Vi € J, (v^, ^-{xa) = 0}- Assume 

that for all i E I, /(xj,x) < 0. XTien /or a// 5 > 7 

lim F(T e K xt A < (a + 5) log if) = 1, 

K—>oo 

where a := sup ie/ (l/|/(xi, x)|). 
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Results as Theorem 4.1 (i) are often called "problem of exit from an attracting domain" 
and can be solved using classical large deviation tools |16| . In the next section, we will 
need a stronger (and non standard) large deviation result for a (small enough) perturbed 
version of the piecewise deterministic process (v^ ,~R K (t),t > 0). Therefore we provide a 



specific proof of the next result based on the Lyapunov functionals used in Theorem 3.2 



Proposition 4.2 Let xi, . . . ,x n be distinct traits and assume that the support of Vq is a 
subset of {x±, . . . , x n }. Consider the population process (y , R^) with perturbed resource 
dynamics: 

dB ^ t ± = 9k -R*-Rf (v«, m ) + ajf (t), (4.1) 

where a£ is a (J-t)t>Q-predictable random process bounded by a constant r\. 

Then for any e > small enough, there exists rjo > 0, V e > and e" < e such that, if 
rj < t]q and (i>q , R^(0)) G B e n(x), the time of exit T ex a of (v]^, R (t)) from B £ (x) is 
bigger than e VeK A T~£ u t with probability converging to 1. 



Proof of Theorem 4.1 By Theorem 3.2 and the continuity of the flow of solutions to 



the chemostat system, for all e" > 0, 

T e ii := sup T e "(u,R) < +oo, 
(u,R)eS 

where T £ « (u, R) is the last entrance time in B £ n /2(x) of the solution of CH(n; x) with initial 



condition (it, R). Then, applying Theorem 2.3 on the time interval [0, T £ »], we obtain that 
P(-Ar-) — > 1 when K —> +oo, where 



Ak :={(^„,R K (T £ »))eiV(x)}. 



Then, choosing e" as in Proposition 4.2 and setting T e = T £ », Point (i) follows from the 
Markov property and Proposition 4.2 (with a K (t) = 0). 

For all i G /, on the event Ak, during all the time interval [T e , e VeK AT^ ut ], the birth rate of 
each individual of trait X{ is bounded from above by rjk(xi)Rk + Ce for some constant C. 
Therefore, the number Z\ of individuals with trait X{ at time T e +t is dominated on the event 
Ak by a continuous time binary branching process Z\ with birth rate r]k(xi)Rk + Cs, 
death rate d{xi) and initial condition Zq = Z % T = K (z/|T , l{ Xi }) < s" K. 

It is well-known (cf. e.g. [TJ p. 109]) that 

-f(xi;x) - Ce 



F(Z l t = | Zl = 1) = 1 



d{xi) exp((-f(xi] x) - Ce)t) - ^ fc Vk{xi)Rk + Ce 



Choosing e small enough for — f(xi; x) > 2Ce and using the fact that ¥(Zl > | Zq = j) = 
1 — F(Zl = | Zq = iy, it is immediate to check that, taking tx = (l/\f(xi; x)| +5)logK, 

lim sup P{Zl K >0\Z G \=j) = 

K->+oc 0<j<Ke 

if e is small enough. Point (ii) then easily follows. □ 
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Proof of Proposition 4.2 In all the proof, we shall denote by u and R the vectors u(x) 
and R(x). 



Let us first observe that, because all the events involved in Proposition 4.2 are T t k 

. . mut 

measurable, it is sufficient to prove Proposition 4.2 for the population process where all 
birth rates with mutations are set to 0. Let us still denote this process by (v^ , (t))t>o- 
Then, at each time t, the set of living traits in the process is a subset of {x%, . . . , x n } and the 
model reduces to a birth and death Markov chain in {-^ r L + ) n for (Nf-(t) := (u^ , t Xi ), 1 < 
i < n )t>o> coupled with (4.1). 

Recall the definition of the Lyapunov functional H (u, R) for the chemostat system in the 

H(u, R) = G(u, R) + jG(u, R), 



proof of Theorem 3.2 



where G is defined in (3.13) and G in (3.17), and 7 > can be arbitrary provided it is 
small enough. 

Set I := {i : Ui = 0} and J := {i : d[xj) = Y^kVn{ x i) R k\- Note that J c C I but in 
general, J c 7^ /. Note that G is linear in the variables (ui,i £ I) and clearly a strictly 
convex function of the variables (ui,i I) and R with diagonal Hessian matrix. Hence, 
since G(u, R) is a quadratic form, making 7 smaller if necessary, 



\u - u\\ 2 + \\R - R\\ 2 < 

i0 



Ui - Ui 



< Ci [H(u,R) -H(u,R)) < Cl(^\ Ui - Ui \ 2 + ^Ui + \\R- R\\ 2 ^, (4.2) 



i0 iel 
for some constant C\ > 1 and for all (u,R) close enough to (u,R). 



By (|3_19j), if (u(t),R(t)) is any solution of CH(n;x), 
d 



dt 



H(u(t),R(t)) < -C 2 \\R-R\\ 2 -2C 2 Y, ^Vk(xi)(ui-Ui) - 2C 2 ^ 



Ui 



k=i \i=i 



< -C 2 \\R-R\\ 2 -C^ (2>fc(*i)( 



Ui - Ui 



fe=i \i<=j 



c 2 y2 u i 



for a positive constant C 2 if (u,R) close enough to (u,R). Due to Assumption (3.2), as 
seen in the end of the proof of Theorem 3.2 the second term of the right hand side is zero 
iff u = u. Therefore, introducing 



C 3 := inf I ^ ( S }2 l r lk{xi){ui - u^ : (w^gj £ s.t. y~] \uj - Uj\ 

k=i ieJ ieJ 



l} >0, 



one deduces that 

d -H(u(t),R{t)) < -C 2 \\R-R\\ 2 -C 2 C 3 Y,\ u i-Ui\ 2 - 



dt 



C 2 J2ui< -d [\\u - u\\ 2 + \\R - R\\ 

(4.3) 
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for some positive constant C4. 

In order to keep notations simple, we shall denote by H(v^ , R K (t)) the function 

H(N?{t),...,N*(t), <(*),..., 
Then, it follows from (|2.9|) that 



J 



1 



{!<«(",-), e<(l- AtK p( a;^ K_)))Efe =1 %(a• ^ K-))^^ fc ( s )} A^ l( ^^s ' d ^ d6, ) 



+ / / / / ( // ( ,f_ + ^i*,l^( a )) - n{y*,B«{B)) 




'0 ^N^O 



i /• /"oo 




'0 JNJO 
ft 



l{i<nK_), 9<m(a i ( I /._),fc)Atirp(a!i(i'.-)) £!=i ^(^K-))^)}^^, d*, d0, ttfl) 

A . - i {i < n( „._), e < diXi(Us _ m N 3 (ds,di,de) 



^ 7o <9i4 

Let T ex it denote the first exit time of {y K (t) , R K (t)) from i? e (x). Assuming i < T ex a, one 
can make a second order expansion of the quantities of the form 

h(v«±^R k (s))-H(v«,R k (s)) 

appearing in the previous equation. Introducing the compensated Poisson point measures 
iVj := Ni — qi for i = 1,2, 3, we obtain, for all t < T ex i t , 



H{yf, R K (t)) = H{v«, R K (0)) + M t K 
n r-t 



+ E / Q-(^,R K (s))N t K ( S )[-d(x l ) + J2vk(x i )R^(s))ds + 
i=l Jo 1 \ k=l ) 

+ E/ § lfk (^,R K (s)){g k -R^(s)(l + (uf,r lk )))ds + 0(r l t), (4.4) 



where is the local martingale 



1 {<<n(i/._), e<(i-MKP(*i(^.-)))ELi%(«i(«'.-))fl»W} Jv i( ds ' <i *' tie ) 

l{i<»(i/._), e<m(* < (i/._),ft)MjrK^(^-))E*=i»7fc(^(^-))«»W}^ 2 ( da ' d *'^' < ^) 



/o vn^o 
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Assuming e small enough, combining (4.2 ), (4.3 ) and (4.4 ), we obtain for all t £ [0, TAT ex [ t ] 
and for K > 1/r) 



N K (t)-u\\ 2 + \\R A (t)-R\\ 2 < C_ 



CA J2\Nf (0)-^| 2 +Vivf (0) + ||i^(0)-i?|| 2 )+ sup \M t 

' te[o,T] 



K 



HI 



iel 



N K (s) - u\\ 2 + \\R K (s) - R\\ 2 - C 5 rj) ds (4.5) 



for some constant C 5 > 0. Therefore, before T A T exi t, \\N K (s) - u\\ 2 + \\R K (s) - R\\ 2 
cannot stay larger than 2C577 on a time interval larger than 



Ci(£i*j " Ui\ 2 + Eiel ^(0) + WHO) - R\\ 2 ) + su Pie[0 ,r] \M< 



tK 



(4.6) 



This implies the following lemma. 



Lemma 4.3 (1) Fix T, 77 > 0, define T v as in (4-6) and let S v denote the first hitting time 
of the 2C^rj-neighborhood of(u,R) by (N K (t),R K (t)). Then, on the event 



we have 



sup (\\N K (t) - uf + \\R K (t) - R\\ 2 ) < dC A C 5 r]T n < 



Tv - TA 2C 1 C 4 C 5 7 ] j' 



->K 



(4.7) 



and < T v A T ex u a.s. 
(2) In addition, on the same event, we also have 

tKu\ „-7l|2 , II r>Ku\ Dl|2 



sup 



N (t) - u\\ 2 + \\R n (t) - R\\ 2 ) < dCiCwpr, + T) < - + ddC^T. 



t<E[0,TAT exU ] 

Therefore, under the additional condition that 



2C1C4C5T 



(4i 



we have T er ji > T. 



We are also going to use exponential moment estimates on the martingale . 

Lemma 4.4 For all a > and T > 0, there exists a constant V a ^ > such that for all 
K large enough 

P( sup \M t K \ > a) < exp(-KV a>T ). 

Proof Since we are only dealing with events occurring before T ex it, we can assume with- 
out loss of generality that the martingale has all its jumps bounded by a constant 
times 1/K, and that all its jump times are (some of the) jump times of Poisson point mea- 
sures which occur at a rate bounded by a constant times K (i.e. the Poisson point measures 
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N^ds^^dO^iKKu^OKrfjg, N 2 (ds, di, d0, dh)li< K u* ,e<rm and N^ds.di^dO)!^^^^, 
where u* := 1 + maxjWj). Therefore, the previous result is a quite standard consequence 
of properties of exponential martingales for pure jump processes. For example, this is a 
consequence of the proof of Proposition 4.1 in |17j . where it is proved that there exists a 
constant C > independent of T and e such that for all a > 

P( sup |Mf| >e) < 2exp(-Kae + CKTT(Ca)), 
Me[o,T] ' 



where t{x) = e x — 1 — x. Since r(x) ~ x /2 when x — >• 0, Lemma 4.4 then follows by 



choosing a > small enough. □ 

Let us now introduce two parameters e' and e" such that e" < e'/2 < e' < e, to be 
determined later. Let To := 0. For all k > 1 such that r^_i < T ex ; t , we define the stopping 
times 

r' k := mf{t > r k ^ : ,K K (t)) ^(x)}, 

r fc := inf{t > r' k : (^f,R X (t)) G fl^x) or (^f,R X (t)) B 6 (x)}. 



Then, combining the strong Markov property with Lemmata |4.3| (1) and |4.4[ it is not 
difficult to prove (cf. that, for convenient choices of e' and e", setting rj = e" /2C§, 
T = 2C\£ 12 /C^C^ri and a = Cie' 2 , these exists a constant V £ := V a> T such that 

sup P(n < T exit ) > 1 - e~ KV * 

(,«R«(0))eB £ ,(x) 

and thus, for all A; > 1, 

sup P(r fc < T exit | r fc _! < T exit ) > 1 - e"^^. 

(^,R^(0))es e ,(x) 

Therefore, if fc ex it denotes the unique integer k such that = T ex ; t , then fc ex it i s larger or 
equal to a geometric random variable of parameter e~ KVs . Hence fc e xit > e^^^ 2 with a 
probability converging to 1. Therefore, in order to complete the proof, it only remains to 
check that for all k, tL — T~k-\ > T' with positive probability for some T" > 0. Thanks to 
the strong Markov property, this is implied by 

inf F(t[ > T') > 0, 

«,R^(0))efJ £ „(x) 



which is a consequence of Lemmata 4.3 (2) and 4.4, replacing e by e'/2 and choosing 
T' = 2C\E n2 jC^C^r\ and a = C\e" 2 (in view of the choice of r\ above, this can impose to 
reduce e" so that 16C 2 e" 2 < e' 2 , but this only changes the constant V £ above), and the 



proof of Proposition 4.2 is completed. □ 



5 Rare mutations and evolutionary time scales 

We now go back to the stochastic model of Section [2] Our goal is to describe the effect of 
the random mutations on the population under a scaling of large population and rare and 
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1 



K = 300, p = 0.0003, a = 0.06, a = 1/4 K = 300, p = 0.0003, a = 0.06, a = 1 

Figure 5.1: Simulations of the individual-based model with rare mutations for two differ- 
ent values of the parameter a. Upper panels: time evolution of the trait density in the 
population. Lower panels: time evolution of the resources concentrations. 



small mutations. We give three convergence results corresponding to three evolutionary 
time scales. These results are similar to the ones obtained the Lotka-Volterra competition 
case (cf . |8] ) and the new ingredients needed in the proofs have been already obtained in the 
previous sections. Therefore this section is devoted to the statement of the results and their 



application to the example introduced in the Section 2.2 and the rigorous modifications of 
the proofs are moved to Appendix. 



5.1 



Convergence to the Polymorphic Evolution Sequence for chemostat 
system 



As shown in Theorem 2.3 mutations have no influence in the limit of large K and small 
\iK on the original time-scale of the population process. We show that evolution proceeds 
on the longer time-scale of mutations ^ * — , for which the population process converges to a 
pure jump process describing successive mutant invasions. After each mutation, disadvan- 
taged traits are driven to extinction due to the competition induced by the deterministic 



chemostat systems. This convergence to a pure jump process is illustrated by Figure 5.1 



The limit process takes values in the set of coupled equilibria of population measures and 
resources, defined by 



Mr 



^^Ui(x)6 Xi , R(x) J ; n > 1, X\, . ■ ■ , x n E X coexist 



vi=l 



where u(x), R(x) were defined in Theorem 3.2 Evolution proceeds by jumps of the support 
of the population measure. 
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We introduce the following assumption: 

For all (xi, . . . , x n ) G T> n , for Lebesgue almost every trait x n +i G X, with the 
notation x = (xi, . . . , x n , x n +i), we assume that for all i G {1, • • • , n + 1}, 

(5.1) 

either Uj(x) > or — <i(xj) + ^J%(xj)i?fc(x) < if Uj(x) = 0. 

fc=i 

In particular, this assumption means that, with the previous notations, defining / = 
{i G {1, . . . ,n + 1} : -d(xi) + Y7k=l r lk{xi)Rk{^-) < 0}, the traits (xj)^/ coexist. 

This assumption implies that we will always be in a clear-cut, non degenerate situation: 
At equilibrium after the invasion of a new mutant trait x n +i in a population with traits 
xi, . . . ,x n , for all of these traits, either it survives or the corresponding growth rate (and 
eigenvalue of the Jacobian) is strictly negative. This allows to use Theorem 4. 1 (ii) . 



, x n G X and 

x " ' ,,,X A - h v itj(x) in probability for all 1 < i < n. Assume 

also that R¥(0) — > flfc(x) in probability when K — > +oo for all 1 < k < r. Assume finally 
that 



Theorem 5.1 Assume (2.61,(2.71,(3.11,(3.21,(5.11. Take coexisting x\ 
assume that v§ = Y17=i u f^Xi with uf 



w>o, iog^:< 



i 



K 



<C exp(VK), as K — > oo. 



(5.2) 



Then, ((vfj KuK ,~R K (tj K /j,f.));t > 0) converges to the A4.Q-valued Markov pure jump process 

{(A1,n a (t));t > 0) de/med as follows: A = E"=i ^( x )^, = R(S?w( A t CT )) w/iere 

Supp(Af) is the support of Af, and the process A a jumps for all j G {1, . . . , n} 



n 

from } Jii(x)5 Xi to 

i=l 



n 

E 

i=l 



ttj(xi, . . . , X n , Xj + /l)^ + U n -fl(xi, . . . , X n , Xj + h)5 X j+h 



with jump measure 



[f(xj + /i;x) 



vfc=l 



ELi%( a; i + /i ) i2 fe( x 



-m a (xj, h)dh. 



(5.3) 



T/ie convergence holds in the sense of finite dimensional distributions on Mf, the set of 
finite positive measures on X equipped with the topology of the total variation norm. 

Following [8j, we call this process Polymorphic Evolution Sequence (PES). The main steps 
of the proof are recalled in Appendix [A| Note that the assumption (5.2) is quite natural 
in view of Theorem |4.1| and ensures that mutation occur after competition eliminates 



disadvantageous traits (Theorem 4.1(h)) and before the population densities drift away 
from equilibrium (Theorem 4.1 -(i)). 

As will appear in the proof, we may interpret the transition rates of the PES as follows: 
starting from an equilibrium population X^=i Uj(x)6 Xj , the process waits an exponential 
time of parameter 



n 

E 



"3i 



\k=l 
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interpreted as a mutation time. The trait of the parent is selected in the population as xj, 
where / is a random variable in { 1 , . . . , n} with distribution 

p/j _ ^ _ P( x i) (ELl W x )) «< 



E"=i P(»i) (ELi Vk(xj)R k (x)) uj ' 
This distribution favors the traits with higher resource consumption. 

Then, the mutant trait is given by xj + H, where H is distributed as m a {xi,h)dh. An 
actual jump occurs if the mutant population, initially composed of a single individual, 
does not go extinct and invades. A comparison argument between the mutant population 
size and branching processes allows to compute the invasion probability as the survival 
probability of a branching process, given by 

[f(x I+ H; X )] + 



After the invasion, the new support of the process At is given by the set of traits with 
nonzero densities in the vector of densities u(x\, . . . , x n , xi + H). 

5.2 The limit of small mutations: the canonical equation 

Until the end of the paper, we assume by simplicity that the trait space is one-dimensional, 
i.e. X C E. 

Our goal is to study the PES under an additional biological assumption of small mutations 
(<7 —> 0), which is standard in this context |19[ [T2| I13| . We prove that when a tends to 
zero, the PES converges on the time scale 4j, to the solution of a (deterministic) ODE, 
called canonical equation of adaptive dynamics, or simply canonical equation. 



By Theorem 5.1 we get Af = u{Xf)5x^ as long as there is no coexistence of two traits in 
the population. Since m r7 (x, h)dh = -m(x, -)dh, the pure jump Markov process (Xf,t > 
0) has the infinitesimal generator 

A a <j)(x)= / (<j)(x + ah) - (p(x))[g(x + ah;x)] + m(x,h)dh, (5.4) 



x 



where 

f{y,x) 



p{x) ^ri k (x)R k {x) u(x) — 
\k=i J 



Vk(y)Rk(x) 



The jump process (X£, t > 0) is often called "trait substitution sequence", or TSS [19J. 
By Proposition |3.5[ the first time of coexistence of two traits in the PES is given by 



r a :='m£{t > : |Supp(Af)| = 2} 

= inf{t > : /(Xf ; X?_) > and f(X?_; X°) > 0}. 



By Corollary B.2 coexistence may only occur in the neighborhood of points x* such that 



dif (%*',%*) = 0. This leads to the following definition. 

Definition 5.2 We call a trait x* £ X an evolutionary singularity if d\f(x*; x*) = 0. 
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Since g(x;x) = 0, the form of the generator (5.4) suggests to scale time as t/a 2 in order 
to obtain a non-trivial limit generator when a — > 0. 



Theorem 5.3 (Theorem 4.4 of [8j) Assume (2.6), (2.7), (3.1), pL2|), (5.1) and that A% 



u(xo)5 Xo where xq is not and evolutionary singularity. Let x(t) be the solution of 
dx{t) 



dt 

such that x(0) = xq. Then, 
(i) For any T > 0, 



h[hdig(x(t);x(t))} + m{x(t),h)dh. (5.5) 



limP(r CT > T/a 2 ) = 1, 



CT->0 



and there exists o~o such that for all a < ao, the process (X^^jt £ [0,T]) is a.s. 
monotone. 

(ii) For any T > 0, the process (A^^it £ |P>^1) converges as a — > to the deterministic 
process (u(x(t))5 x ^,t £ [0, T]) for the Skorohod topology on D([0,T],Mq), where 
A4q is equipped with the weak topology. 



Equation (5.5) is known as the canonical equation of adaptive dynamics [12J. 

This result follows from the convergence of the TSS process X°, 2 to x(t) when a — > 
Thm. 4.1], from the fact that mutation jumps in the trait space are bounded by aDiam(X) 
(by the definition of the jump measure m a ) which converges to when a — > and from a 



careful study of the solutions of Equation (5.5). 



Remark 5.4 In the case when m(x, ■) is a symmetrical measure on K for all x £ X , 

heuristically introducet 

K(x(t))d ig (x(t);x(t)), 



Equation (5.5) gets the classical form, heuristically introduced in [12j, 

dx(t) 1 



dt 2 

where K{x) is the variance of m(x,h)dh. 

5.3 Small mutations and evolutionary branching 



We will now introduce the last tools to predict the different behaviors observed in Figure 2.2 
and in particular to characterize the diversification phenomenon of evolutionary branching. 
Since it is not captured by the canonical equation, evolutionary branching can only occur 
on a longer time scale and in the neighborhood of equilibria of the canonical equation. Let 
x* := lim^+oo x(t) ( well-defined since X is compact). We make the additional assumption 



x* is in the interior of X, / m(x*,h)dh > and / m(x*,h)dh > 0. (5-6) 

JR- Jr + 

This assumption means that mutations are always possible from x* . Then in view of the 



canonical equation (5.5), d\f(x*;x*) = 0, i.e. x* is an evolutionary singularity. 
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The linear stability condition for the equilibrium x* implies that 

d 11 f(x*;x*) + d 12 f(x*;x*)<0. 
Differentiating twice the equation f(x;x) = implies that 

dnf(x; x) + 2d 12 f(x; x) + d 22 f(x; x) = 0, Vx e X , 

and so 

dnf(x*;x*)<d 22 f(x*;x*). 

We shall leave the degenerate case d 22 f(x*;x*) = duf(x*;x*) for further studies, and 
assume below that 

d n f(x*;x*) <d 22 f(x*;x*). (5.7) 

Let us recall the definition of evolutionary branching introduced in [H]. 

Definition 5.5 Fix a > and x* an evolutionary singularity. For all rj > 0, we say that 
there is r]-branching at x* for the PES A a if 

• there exists t\ > such that the support of the PES at time t\ is composed of a single 
point belonging to [x* — r], x* + rf\; 

• there exists t 2 > t\ such that the support of the PES at time t 2 is composed of exactly 
2 points separated by a distance of more than r//2; 

• between t\ and t 2 , the support of the PES is always a subset of [x* — r], x* + rj], and 
is always composed of at most 2 traits, and has nondecreasing (in time) diameter. 

This definition only considers binary evolutionary branching. The next theorem proves 
that, in the neighborhood of an evolutionary singularity, evolutionary branching of a mo- 
nomorphic population into a n-morphic population with n > 3 is impossible (at least when 
the trait space is one-dimensional) . Note that the notion of evolutionary branching requires 
the coexistence of two traits, but also that these two traits diverge from one another. 



Theorem 5.6 (Evolutionary branching criterion) Assume (2.61,(2.71,(3.11,(3.21,(5.1) 



and (5.6). Assume also that Aq = u(xq)5 Xo and that the canonical equation 



5.5) with ini- 



tial condition xq converges to an evolutionary singularity x* in the interior of X . Assume 



finally that x* satisfies (5.7) and 

d 22 f(x*;x*)+d u f(x*;x*)^0. (5.8) 
Then, for all sufficiently small rj, there exists ctq > such that for all a < o~q, 

(a) if duf(x*\ x*) > 0, F(rj- branching at x* for K a ) = 1. 

(b) if d\xf(x*\ x*) < 0, ¥(rj -branching at x* for A a ) = 0. Moreover, 

F(Vi > 0°{x*), Card(Supp(An) < 2 and Supp(Af) C (x* - r ) , x* + rj) ) = 1, 
where 

9°(x*) = mi{t > 0, Supp(A£) n (x* - r), x* + rj) ^ 0}. 
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This criterion appeared for the first time in |19[ Section 3.2.5] with an heuristic justification. 
Locally around x*, one of the two following events can occur almost surely: either there is 
binary evolutionary branching or not. Coexistence can occur in both cases, but our proof 
shows that, after coexistence, the diameter of the support of the PES is a.s. non-decreasing 
(resp. non-increasing) in the first (resp. second) case. 

Remark 5.7 It can be proved using Proposition \B. 3| (i) that 

• if d\\f{x*\ x*) +d22f(x*;x*) > 0, then for all neighborhood U of x* in X, there exist 
x,y £U that coexist and 

• if d±if(x* ; x*) + d22f(x* ; x*) < 0, then there exists a neighborhood U of x* in X such 
that any x,y 6W do not coexist. 



Comparing this result with the criterion of Theorem 5.6, one has of course that evolutionary 



branching implies coexistence, but the converse is not true. 



As in [8], the proof of Theorem 5.6 is based on local expansions of the fitness functions 



f(y;x) and f(z;x,y) in the neighborhood of x*, given in our case in Propositions B.l 
and |B.3| in the Appendix. The main steps of the proof of Theorem |5.6| are given in 
Appendix [C] 

5.4 Back to our example 



Let us come back to the example developed in Subsection 2.2 The death rate of an 
individual with trait x has the form d(x) = \ + ax 2 with a > 0. Let us compute in this 
case the quantities we are interested in. By symmetry considerations, we have that x* = 



is an evolutionary singularity. It then follows from (3.20) that u(0) = 3 and from (3.9) 
that 

dnf{P,0)--d (0)+ 1+??i(0) - (0) -l"2a. 
For symmetry reasons again, we also have u'(0) = 0, from which we deduce that 

<9 2 2/(0;0) = —■ 

Differentiating twice the relation 

+ ax 2 J (l+{x-l) 2 u{x)){l+(x+l) 2 u{x)) = (x-l) 2 (l+(x+l) 2 u(x))+{x+l) 2 (l+(x-l) 2 u 



we obtain u"(0) = -4 - 16a, which give <9 2 2/(0; 0) = 4 + 2a. 

Therefore, the inequality ^22/(0; 0) > f?n/(0; 0) is always satisfied, and there is evolution- 
ary branching (resp. no evolutionary branching) at x* = if a < 1/2 (resp. a > 1/2). In the 
case where a < 1/2, the death rate does not increase too much as the trait x deviates from 
0. So the population is better off dividing into two specialized populations, each consuming 
better either resource 1 or resource 2. In addition, since 9n/(0;0) + 02 2 /(0;0) = 5 > 0, 
coexistence is always possible in the neighborhood of x* = 0. Therefore, the simulations 



of Figures 2.2 and |5.1|are consistent with Theorem 5.6 and Proposition 5.7 
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Appendix 

A Sketch of the proof of Theorem |5.1 



We now give the general idea of the proof of Theorem 5.1 extending the biological heuristics 



of |19j . rigorously proved in [5] in the case of Lotka-Volterra competition. 

Let us roughly describe the successive steps of mutation, invasion and competition. The 
two steps of the invasion of a mutant in a given population are firstly the stabilization 
of the resident population before the mutation and secondly the invasion of the mutant 
population after the mutation. 

Fix t] > 0. In the first step, assuming that n traits x\, . . . ,x n coexist, Assumption K ^ <C 
e VK and Theorem 4.1 ensure that the population densities ((vf, l{ Xl y), . . . , (f/ 4 ", l{ Xn y)) 
and the vector of resources R K (t) belong to the //-neighborhood of (u(x),R(x)) with 
high probability for large K until the next mutant y appears. Therefore, until T^ u ^^ the 
population densities are roughly constant and the rate of mutation is approximated by 
^KP{xi) Y2kVk(xi)Rk( x )Kui(x). More precisely, the next lemma can be proved with the 
same arguments as in Lemma 2 (b) and (c) in [5]. 

Lemma A.l Let Supp(uq) = {xi, . . . ,x n } that coexist. There exists Eq > such that, 
if (u^, RA(0)) G B £o (x), then, for any e < e , 

lim Pfc > log if; Vt G [log K, T%J , {yf ,R K (t)) G S e (x)) = 1, 

n r 

K ^KT* ut => Exp(y2p(x j )y2ri k (x j )R k (x] 

A^-OO \ * — * *• — rf 

3=1 k=l 

and 

lim F(at time '/',(),,,. the mutant is born from trait .r, ) - ' —/—J. 



"""" E"=iP(^?)£Ll%( x j)^fc( x ) 

for alii G {1, . . . , n}, where =^ stands for the convergence in distribution of real random 
variables and Exp(a) denotes the exponential distribution with parameter a. 

In the second step, we divide the invasion of a given mutant trait y into 2 phases shown 
in Fig. |A.1| 

The first phase stops either when the mutant population gets extinct or reaches a fixed 



small density rj > (at time ti in Fig. A.l). During all this phase, the mutant density is 



small, and so, using the perturbed version of the large deviation result (Proposition 4.2), 



we can prove that the resident population stays close to its equilibrium density u(x) and 
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population size 



u(x) 



V 



t 2 



Figure A.l: The two steps of the invasion of a mutant trait y in a monomorphic population with 
trait x. 



the resource concentrations stay close to R(x). Therefore, similarly as in the proof of The- 



orem 



4.1 ■■(ii), the number of mutant individuals can be compared with branching processes 



with birth rate r/k(y)Rk(x-)±Ce and death rate d(y). The growth rate of this branching 
process is close to the fitness f(y; x), which hence describes the ability of the initially rare 
mutant trait y to invade the equilibrium resident population with traits x\, . . . , x n . If this 
fitness is positive (i.e. if the branching processes are super-critical), the probability that 
the mutant population reaches density rj > at some time ii is close to the probability 
that the branching process reaches r]K, which is itself close to its survival probability, i.e. 
close to [/(y;x)] + / (^^Vkiu)^^))) when K is large. In addition, the comparison with 
branching processes shows that this first phase ends before a time of order log K similarly 
as in the proof of Theorem 4.1 - (ii) (see the proof of Lemma 3 in [5]), and so the assumption 



log if < 



Kuk 



(A.l) 



ensures that no new mutation occurs during this first phase with high probability. 



If the mutant population invades (i.e. reaches the density 77), the second phase stops 
when the population densities and the resource concentrations reach B„(x), where x = 
(x±, . . . , x n , y), and when the traits i such that Uj(x) = go extinct in the population (at 
time ti in Fig. |A.l ). Theorem 4.1 (i) and-(ii) ensure that this phase is completed with 
a probability converging to one after a time of order log if. Note that Assumption (5.1) 



is required to be able to apply Theorem 4.1 (i) for almost all mutant trait y born in 
a population with traits x±, . . . ,x n . Again, the assumption (A.l) ensures that no new 



mutation occurs during this second phase with high probability. 

These two phases are summarized in the following lemma, which can be proved as Lemma A. 4 
in 0. 
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Lemma A.2 Let Suppiy^) = {x\, . . . , x n , y} where x\, . . . , x n coexist and y is a mutant 



trait that satisfy Assumption (5.1). We shall denote x n +i = y and x = (x\, . . . , x n , x n +\) 
for convenience. We define 

rf = inf{t > : {vf , K K "(t)) G £ 6 (x) and Vi a.t. 14 (x) = 0, (z/f , = 0} 

rf = inf{t > : (i/f , l {y} ) = and (z^f , R5(i)) G 5 e (x)}. 

Assume that (v^,ls y y) = 1/K (a single initial mutant). Then, there exists Eq such that 
for all e < e and if (z/jf , R^(0)) G B e (x), 

lim P(rf <7 f) = -Jfe^!± Hm P(r K <T K ) = 1 [/(*;*)]+ 



and Vn>0, lim P ( if A r 2 ^ < — — A T^. ,] = 1. 

Combining all the previous results, we can prove as in [5] that for all e > 0, t > and 
r C X measurable, 

lim F(A £n (t, T)) = pfSupp(A^) C T and has n elements') (A.2) 

where the event 

^4 £jn (t,r) := ^Supp(^ K ^ K ) C r has n elements that coexist, say x\, . . . ,x n , 

and VI < i < n, \(v$j Kllk , l{ Xi} ) - u;(x)| < e} 



and (A%,t > 0) is the PES defined in Theorem 5.1 
The proof ends as in [5] . 

B About the sign of the fitness functions 



In the PES of Theorem 5.1 the success of a mutant invasion is governed by the sign of its 



fitness. Our goal in this section is to study the fitness of mutant traits in the neighborhood 



of coexisting traits. These results are used in Section 5.2 to study the local direction of 



evolution in the PES, and the phenomenon of evolutionary branching. 

We will assume in all what follows that the traits are one-dimensional, i.e. X C M. 

Proposition B.l For all x* = (x*,...,x*) G T> n , we have for all 1 < % < n, when 
x = (x\, . . . , x n ) G V n — > x* and y — > x*, 



/(y;x) = (y - x^ 



no + y ° kV ' k{ f\ , . + (i) ) . 
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Proof Using the relation /(xjjx) = 0, we have 

9k(vk(y) - Vkjx i)) 



ft \ 7/ \ I if \ | ^ ^ k\y - i k^i 

= -(2/ - x.) ^ + (y - + (y - **) £ H^E^feS^- 



Proposition B.l then easily follows from Lemma 3.6 □ 



If we think of y as a mutant trait born from the resident trait x* in the resident population 
of traits x* = x, this result gives the sign of the fitness function when mutations are small, 
i.e. when \y — x*\ small, except when the first-order term is zero. Note that this first-order 
term is given by the derivative of the fitness function with respect to the first variable at 

(x*,x*)\ 

Observe also that, for all 1 < i < n, /(a^jx) = for all xi,...,x n which coexist. Hence, 
if (xi, ...,x n ) is in V n , we have dif(xi]x.) = -d i+ if(xi;x). 

The differentiability of the fitness function at a point (y, x\, . . . , x n ) G X X T> n is ensured by 
Lemma 3.6 The following result shows that coexistence of n + 1 traits (x, y) with x close 
to x* G D n and y close to one of the coordinates x* of x* is only possible if d\f(x*; x*) = 0. 

Corollary B.2 If x* G D n and dif(x*;x*) ^ 0, then there exists a neighborhood V of 
(x*,x*) = (x\, . . . such that V PlV n+1 = 0. 

Proof Assume that for all neighborhood V of (x*,x*), V n £> n +i 7^ 0- Since P n is 
an open subset of X n , there exists x close enough to x* and y close enough to x* such 



that x G V n , (xi, . . . , Xi-i,y, x i+ i, ... ,x n ) G £>„, (x, y) G P n+i and, by Proposition B.l 



f(y; x) > and f(xf,Xi, . . . , y, a?i+i, . . . , x n ) < 0. By Theorem 3.2 , this means that 
u(x, y) = (u(x),0), which contradicts the fact that (x, y) G "D n +l- ^ 

In the neighborhood of points x* G T) n such that d\f(x*,x*) = 0, we have more precise 
results on the sign of the fitness. 

Proposition B.3 Fix x* = (x*, . . . , x* n ) G V n and 1 < i < n such that 

(i) When (x\, . . . ,x n ) — )• x* and y — > x*, 

/(y; x) = i(y - ^) f (a + o(l))(y - x*) + + o(l))( Xj - x*) 

where 

a = -d»( X *) + Y " k4( f\ , 6, = a + 2 V E ^ ^iM^M^ 

^1 + Er=i%(4)^(x*)' ^ [1 + Er=i%(4)^(x*)] 2 
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and for j 7^ i, 



2 sp 9kV k (x* ) Ya=i d xj (Vk {xj)ui (x) ) (x* 

[i+Er=i%(4)^(x*)] 2 



fc=i 



(ii) Assume that for all neighborhood V of (x*, x*) in X n+1 , V n £> n +i 7^ 0. Then, when 

x — > x* and y — )• x* such that (x, y) £ T) n+ i, 

Ui (x, y) + n n+ i (x, y ) — > Uj (x* ) . 

(iii) Assume that for all neighborhood V of (x*, x*) in X n+1 , V fl fn+i 7^ 0. Then, when 

x ->■ x* ; y ->■ x* and z ->■ x*, 



f(z;*,y) = -{z - Xi)(z - y)(a + o(l)). 
If one thinks of y as a mutant born from trait Xj in a coexisting population of traits x, 



Point (i) allows to characterize the cases where (x, y) G thanks to Proposition 3.7 



Point (ii) shows that, in this case, the sum of the densities of traits y and X\ is close to 
«j(x*). Point (iii) shows that if a second mutant z is born from either y or Xj, then its 
fitness is positive or negative, depending on the sign of a and on the position of z with 
respect to Xi and y. If z is a mutant trait born from another resident trait Xj with j 7^ i, 
then the sign of the fitness function can be determined using Proposition B.l and Point (i) 
of Proposition B.3| 

Proof Since dif(x*;x*) = 0, we have 



(B.l) 



and therefore 



f(y, x ) = ~{y - x i) / [d'(xi + (y - Xi)u) - d'(x*)]du 



. / ^5Jo[#'+ (y-Xi)n) -nj,(x*)]dn 



1 + T,jVk(x j )u j (x) 



+ {y-x 



Xj)Uj (x) 1 + £\ Vk {x*)uj (x* ) 



(y - Xj)(d (x*) + o(l)) / (xi + (y - x;)u - x*)du 

Jo 

9k(Vki x *) + oi 1 )) fo( x i + (y~ Xi)u - x*)du 



+ {y- xi) 

k=i 



+ (y-xi)J2j 1 



9kV' k (x*) 



^a + E 3 m(x*)u,(x*)f 



(x* -x)* ^V(%(x i )%(x))(x*) + (1) , 



where the differentiability of x 1— y n(x) comes from Lemma 3.6 Point (i) easily follows. 
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Point (ii) can be proved as follows. Recall that, by Lemma 3.6 u is a bounded function on 
X n+l . Let u* be any accumulation point Uj(x, y) + u n+1 (x, y) when x — > x* and y — > x* 
such that (x, y) £ T> n +\- After extracting a subsequence, we may also assume that for all 
j h Uj(x,y) -> u*. 

Passing to the limit in (3.24), we get for all j G {1, . . . , n}, 



Vk{x*)g k 



m=l Vkv^m) ^"m 



d(x*). 



By Assumption (3.2), the unique solution of this system of equations is u(x*) and the 
convergence follows. 

Point (hi) can be proved with a similar computation as for Proposition B.l using (3.24) 
and the fact that f(xf, x, y) = 0, we have 

/O; x, y) = -(z-xi) [ d' (xi+(z-Xi)u)du+(z-Xi) V - 9k ^ + — x i) u ) du 

Jo ~ ^ + 22j=iVk(xj)u j x,y + r] k {y)u n+1 {x,y) 



Then, using the relation 

Z X{ 







y-xi 



f(y;x,y) = -(z-Xi) / d '(xi + (y - Xi)u)du 



+ (z- x-) V gk Jo v k (xj + (y - xj)u)du 

(r[ 1 + YJj=i r lk{xj)u j x,y + rj k (y)u n+1 (x,y)' 



we obtain 



f(z; x, y) = -(z - x^ / [d'(xi + (z - Xi)u) - d'(xi + (y - Xi)u)]du 



+ (z- x-) V 9k k^M^l + ( Z ~ ~ r M^i + (y~ x i) u )] du 
l + Y^ r LiVk(xj)ujX,y + i]k(y)u n+1 (x,y) 



k=l 



= -{z - Xi)(z - y) 
+ (z — X 

and the result follows from Point (ii 



u 

JO 



d"(xi + (y — Xi)u + (z — y)uv)dv du 



w _ v fffc Jq 1 u 77^ (xj + {y - Xj)u + {z- y)uv)dv du 
fr[ 1 + T l ]=iVk(x j )u j x,y + r] k (y)u n+1 (x,y) 



□ 



C Sketch of the proof of Theorem 5.6 



We recall the main steps of the proof of the branching criterion in [H]. 

1. Before the first coexistence time, the support of the PES, given by the TSS, reaches 
(x* — r/,x* + 77) in finite time almost surely (at time 9Z{x*)). 

2. After time 9Z(x*), the (support of the) PES cannot exit (x* — 77, x* + 77) while being 
monomorphic. 
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If d\if(x*;x*) < 0, either there is never coexistence, and then the result is proved, 
or there is coexistence of two traits after a finite time, and then, if two traits x < y 
(say) in (x* — 77, x* + 77) coexist, 



(a) by Proposition B.3 (iii), the only mutants traits z that can invade the dimorphic 
population of traits x and y are such that z G (x,y); 

(b) the three traits x, y and z cannot coexist, and at least one of the resident traits 
x and y goes extinct, i.e. u\(x,y,z) = or U2(x,y,z) = 0. 

This shows that, after any coexistence time, the distance between the two branches 
can only decrease, until one of the branches goes extinct. 

If dnf(x*;x*)>0, 

(a) there is almost surely coexistence of two traits in finite time; 



(b) by Proposition B.3 (iii), if two traits x < y (say) in (x* — 77, x* + 77) coexist, the 
only mutant traits z that can invade the population of traits x and y are such 
that z G" [x, y]; 

(c) once such a mutant trait invades, x, y and z cannot coexist, and the intermediate 
trait only goes extinct, i.e. u(x, y, z) G W + x {0} xW + if x < y < z, or u(x, y, z) G 
{0} x (R* + ) 2 if z < x < y. 

This shows that the two branches survive until one of them exits of {x* — 77, x* +77), 
and the last point of the proof consists in proving that the distance between the two 
branches becomes bigger than 77/2 in finite time almost surely. 

The only steps that require a different proof from Theorem 4.9 of [8 J are Steps 3.(b) 
and 4.(c). Indeed, in [8], these steps were proved using general results on the long time 
behavior of 3-dimensional competitive Lotka-Volterra systems [26] • Here, we do not have 



such general results, but we have Theorem 3.2 



Proof of Step 3.(b) Let us assume that dnf(x*; x*) < 0, and let x < y G (x* — 77, x* +77) 



coexist, i.e. by Corollary B.2 f(x;y) > and f(y;x) > 0. Let also z be a mutant trait in 
(x* — rj, x* + 77) that can invade the resident population of traits x and y, i.e. z G (x, y) 
and f(z; x,y) > by Proposition B.3| (iii). 



Using the relation f(x; y) - f(y; y) = f(x; y) = (x - y) f Q d\f{y + u(x - y))du, we have 
° ( Hx:u) ^ f ud u f( y + u(x-y);y)du. (C.l) 



dx \ y — x 







Since d\\f{x*;x*) < 0, this is positive if x and y are sufficiently close to x* . Hence, since 
x < z < y, we have f(z; y)/{y — z) > f(x; y)/(y — x) > 0, and thus f(z; y) > 0. Similarly, 



f(z;y) > 0. In addition, by Proposition B.3 (iii), if x and z coexist (i.e. if f(x;z) > 0), 



then f(y; x,z) < and if y and z coexist (i.e. if f(y; z) > 0), then f(x; y, z) < 0. Therefore, 
by the characterization of u(x, y, z) of Theorem |3.2[ if f(x; z) > 0, then 



u(x,y,z) = (ui(x,z),0,u 2 (x,z) 
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and if f(y; z) > 0, then 

u(x, y, z) = (o, ui(y, z),u 2 (y, z)j . 

Conversely, in the case when f(x;z) < and f{y;z) < 0, then by Theorem 3.2 again 
u(x,y,z) = (0,0, u(z)). 

Finally, we see that x, y and z never coexist, and that at least one of the traits x or y 
goes extinct in u(x, y, z). Therefore, once two traits in (x* — rj,x* + rj) coexist in the PES, 
then the next jump in the PES reaches a position where either only one trait survives in 
(x* — r/,x* + rj), or two traits survive in (x* — T), x* + rj), closer to each other than before 
the jump. □ 



Proof of Step 4.(c) Let us assume that dnf(x*; x*) > 0, and let x < y G (x* — rj, x* +rj) 



coexist, i.e. by Corollary B.2 f(x;y) > and f(y;x) > 0. Let also z be a mutant trait in 
(x* — rj, x* + rj) that can invade the resident population of traits x and y, i.e. z (x, y) 
and f(z; x,y) > by Proposition B.3| (iii). Let us assume for example that z < x < y. 



Then, it follows from (C.l ) that f(z; y) > 0. Similarly, since by assumption (5.7) d22f( x *; x* 
3n/(x*;x*) >0, 



d_ 

Ox 



ud 2 2f{y; y + u(x- y))du < 



for all x and y close enough to x*. Therefore, f(y;x) > implies that f(y;z) > 0. Hence 



y and z coexist and by Proposition B.3 (iii) f(x;y,z) < 0. All these conditions imply by 
Theorem 13.21 that 

u(x,y,z) = (0,ui(y,z),u 2 (y,z) 



Similarly, if x < y < z, then 



u(x,y,z) 



ui(x, z), 0, u 2 (x, z) 



In all cases, x, y and z cannot coexist, and once two traits in (x* —rj,x*+ rj) coexist in the 
PES, then the next jump in the PES reaches a position where two traits survive, farther 
from each other than before the jump. □ 



> 
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